Dynamical Instability of a Doubly Quantized Vortex in a Bose-Einstein condensate 
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Doubly quantized vortices were topologically imprinted in \F = 1) 23 Na condensates, and their 
time evolution was observed using a tomographic imaging technique. The decay into two singly 
quantized vortices was characterized and attributed to dynamical instability. The time scale of the 
splitting process was found to be longer at higher atom density. 
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Quantum fluids, like superfluid He, electrons in a su- 
perconductor or a Bose-Einstein condensate of atoms, are 
described by a macroscopic wavefunction. This requires 
the flow field to be irrotational, and gives rise to super- 
fluidity and quantized circulation 1]. Atoms in a Bose- 
Einstein condensate, for example, can only circulate with 
angular momentum equal to integer multiple of Ti, in the 
form of a quantized vortex |2|. 

Vortices are excited states of motion and therefore en- 
ergetically unstable towards relaxation into the motional 
ground state, where the condensate is at rest. How- 
ever, quantization constrains the decay: a vortex in Bose- 
Einstein condensates cannot simply fade away or disap- 
pear, it is only allowed to move out of the condensate 
or annihilate with another vortex of opposite circulation. 
Vortex decay and metastability, due to inhibition of de- 
cay, have been a central issue in the study of superfluid- 
ity 0, m ls M S m • I n almost pure Bose-Einstein conden- 
sates, vortices with lifetimes up to tens of seconds have 
been observed 0, EH3 • 

Giving a Bose-Einstein condensate angular momentum 
per particle larger than % can result in one multiply- 
quantized vortex with large circulation or, alternatively, 
in many singly-quantized vortices each with angular mo- 
mentum h. The kinetic energy of atoms circulating 
around the vortex is proportional to the square of the 
angular momentum; therefore the kinetic energy associ- 
ated with the presence of a multiply-quantized vortex is 
larger than the kinetic energy of a collection of singly- 
quantized vortices carrying the same angular momen- 
tum. A multiply-quantized vortex can decay coherently 
by splitting into singly-quantized vortices and transfer- 
ring the kinetic energy to coherent excitation modes, a 
phenomenon called dynamical instability which is driven 
by atomic interactions 0, EJ Q > and not caused by 
dissipation in an external bath. Observations of arrays of 
sing ly-q uant ized vortices in rapidly rotating condensates 
lllj indirectly suggests that the dynamical instabil- 
ity leads to fast decay of multiply-quantized vortices. 
However, the existence of stable multiply-quantized vor- 



tices in trapped Bose-Einstein condensates has been pre- 
dicted with a localized pinning potential or in a 
quartic potential [l5|- Stable doubly-quantized vortices 
were observed in superconductors in presence of pinning 
forces and in superfluid 3 He- A which has a multicom- 
ponent order parameter Recently, formation of a 

multiply-quantized vortex in a Bose-Einstein condensate 
has been demonstrated using topological phases 0, ^| , 
and surprisingly long lifetime of a "giant" vortex core 
has been reported [2fl|. The study of topological excita- 
tion and their stability is an active frontier in the field of 
quantum degenerate gases 0, [2^ . 

In this Letter, we study the time evolution of a doubly- 
quantized vortex state in a Bose-Einstein condensate, 
and directly confirm its dynamical instability by observ- 
ing that a doubly-quantized vortex core splits into two 
singly-quantized vortex cores. The characteristic time 
scale of the splitting process was determined as a func- 
tion of atom density and was longer at higher atomic 
density. 

Bose-Einstein condensates containing over 10 T 23 Na 
atoms were created in the |F = 1, ra^ = — 1) state, trans- 
ferred into an auxiliary chamber |23|, and loaded into a 
Ioffe-Pritchard magnetic t rap generated by a microfab- 
ricated atom chip M [H [m. The wire pattern on the 
atom chip is shown in Fig. da). In our previous work, 
we used a Z-shaped wire trap where changing the sign 
of the axial magnetic field curvature was technically im- 
possible so that we could not trap condensates after im- 
printing a vortex. To overcome this technical difficulty, 
we designed our new chip with separate end-cap wires, 
allowing independent control of the axial magnetic field. 
Typical wire currents were Ic = 1.53 A in the center wire 
and II = Ir = 0.1 A in the end-cap wires, and the ex- 
ternal magnetic field was B z — 450 mG and B x = 5.3 G, 
resulting in a radial (axial) trap frequency f r = 220 Hz 
(fz — 3 Hz) and a distance of the trap from the chip 
surface d = 600 fim. After holding condensates for 2 s 
to damp excitations which might have been caused by 
the loading process, condensates contained over 1.5 x 10 6 




FIG. 1: (a) Wire pattern on the atom chip. A magnetic 
trap is formed by a current Ic flowing through the center 
wire in conjunction with a external uniform magnetic field 
B x . The axial confinement along z direction is generated 
by currents II and Ir in the end-cap wires. Each current is 
controlled independently. A 2 /im thick Au film was deposited 
by evaporation on a thermally oxidized Si substrate and wires 
were patterned by photolithography and wet etching. The 
width of the center wire and the end-cap wires were 50 /im 
and 100 /im, respectively, (b) Imprinting of a vortex in a Bose- 
Einstein condensate. By inverting the z direction magnetic 
field B z , a doubly quantized vortex was imprinted in \F = 1) 
condensates, using topological phases as in Ref. 19]. The 
direction of II and Ir were also reversed to maintain the axial 
confinement. The dashed lines indicate the selective probing 
region for tomographic imaging as described in the text. 



atoms and the lifetime of condensates was ~ 8 s with a 
radio- frequency (rf) shield [27^ . 

Doubly-quantized vortices were topologically im- 
printed in condensates by inverting the axial magnetic 
field, B Zl as demonstrated in Ref. 19]. B z was ramped 
linearly from 450 mG to —460 mG in 12 ms. As B z 
passed zero, the sign of axial field curvature was changed 
by reversing the directions of II and Ir in 1 ms. The 
trap center position and the axial trap frequency of the 
inverted wire trap were matched to those of the origi- 
nal wire trap by adjusting the final values for II and 
Ir. Losses due to nonadiabatic spin flips as B z passed 
through zero reduced the number of atoms in the conden- 
sate after imprinting to about ~ 1 x 10 6 , giving a typical 
healing length £ = 0.4 jam. The lifetime of condensates 
after imprinting was less than 2 s. 

The vortex imprinting process was accompanied by a 
sudden mechanical squeeze in the radial direction and a 
kick in the vertical direction. The radial trap frequency 
is proportional to the square root of the bias magnetic 
field (f r oc l^l -1 / 2 ) and became temporarily higher dur- 
ing field inversion. Additionally, the vertical position of 
the trap center changed as the gravitational sag (ex f~ 2 ) 
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FIG. 2: Decay of a doubly quantized vortex. Axial absorption 
images of condensates after 15 ms of ballistic expansion with a 
variable hold time after imprinting a doubly quantized vortex. 
A doubly quantized vortex decayed into two singly quantized 
vortices. For this data, the interaction strength was an z ~ 7.5 
(see text for definition). The field of view in each image is 
320 /xm x 320 /im. 



changed from 5.1 /im to zero. The Thomas- Fermi radius 
of condensates in the loading phase was ~ 5 /im. After 
imprinting a vortex, the amplitude of quadruple oscilla- 
tion in the axial direction was ~ 20 % of the axial length 
of condensates (~ 600 /im), but there was no detectable 
dipole oscillation in the vertical direction. 

The decay of a doubly-quantized vortex state was stud- 
ied by taking an absorption images along the imprinted 
vortex line after releasing the condensate and letting it 
expand for 15 ms. When we took an integrated absorp- 
tion image, the visibility of a vortex core completely van- 
ished within 30 ms. To reduce blurring due to possible 
bending of the vortex line l2fll . we employed a tomo- 
graphic imaging technique |21j. A 30 /im thick central 
slice of the condensate (see Fig. [Hb)) was selectively 
pumped into the F = 2 hyperfine level with a sheet 
of laser light perpendicular to the condensate long axis; 
the radial profile of the condensate in the selected region 
was then imaged with a light pulse resonant with the 
F = 2 — » F' = 3 cycling transition. In our absorption 
images, the size of a doubly-quantized vortex core was 
typically ~ 40 jam. This tomographic imaging technique 
was crucial for observing the time evolution of vortex 
cores beyond 30 ms. 

A series of absorption images of the splitting process 
is provided in Fig. [2j Images taken just after imprinting 
show a doubly-quantized vortex core of high visibility; 
the visibility of the core decreased with time, an effect 
we attribute to bending of the vortex line [2^ and other 
excitations created during the imprinting process. Later 
in the evolution, the central core deformed into an ellip- 
tical shape and split into two closely-spaced cores. Once 
the two cores were separated by their diameter, they ap- 
peared well resolved in our images. The angular position 
of the two cores was random for each experimental real- 
ization with the same evolution time, so the precession 
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frequency of two cores could not be determined with our 
destructive image technique. 

To investigate the dependence of the instability on the 
mean field atomic interaction, we measured the charac- 
teristic time scale of splitting of a doubly-quantized vor- 
tex core as a function of the atom density. Atom density 
was controlled by removing a variable number of atoms 
with rf evaporation before imprinting a vortex. Images 
were classified as follows: images where the two cores 
were separated by more than one core diameter were 
labelled as "two visible cores"; images with a clearly- 
defined circular central core were labelled as "one core" ; 
images in the intermediate range, where the central core 
was elliptical but the two cores were not resolved, or with 
a bad visibility were labelled as "undetermined" . For ex- 
ample, the images at 62 ms and 75 ms in Fig. and 
Fig.OJa) were classified as "two visible cores", and 50 ms 
in Fig. [2j and Fig. Ufa) an d (c) as "undetermined" . 

Experimental results are provided in Fig. 03 as a func- 
tion of the linear atom density n z (along the conden- 
sate long axis) multiplied by the s-wave scattering length 
a. The rescaled quantity, an z = a J \tlj(r)\ 2 dxdy corre- 
sponds for a cylindrical condensate to the strength of the 
mean field interaction, with ip(r) being the condensate 
wavefunction. Results in Fig. 03 clearly demonstrate that 
a doubly-quantized vortex core splits more slowly as the 
density becomes higher. 

Once the doubly-quantized vortex core split into two 
cores, the distance between the two cores was almost con- 
stant (~ 50 jam) during the further evolution, as shown 
in Fig. Efc). This is evidence that the separation pro- 
cess was driven mainly by the dynamical instability, and 
not by dissipation, which would gradually increase the 
separation of the two cores. Dissipative processes were 
minimized by performing the experiments at the lowest 
possible temperature. Condensates did not have any dis- 
cernible thermal atoms even after extended hold time. 
Furthermore, the energy released by the dissociation of 
the doubly-quantized vortex was ~ 5 nK negligible to 
the critical temperature ~ 240 nK. For the upper bound 
to the temperature of < 100 nK, Ref. 30] predicts that 
dissipative decay time to be « 1.5 s for a single vortex, a 
time scale much longer than what we observed. 

Multiply-quantized vortices in a harmonic potential 
are predicted to spontaneously decay into other states 
even in the absence of dissipation and external pertur- 
bations p. In the Bogoliubov framework, which is be- 
lieved to well describe quantized vortices in one compo- 
nent condensates, the dynamical instability manifests as 
the existence of excitation modes with a complex eigen- 
frequency. The nonvanishing imaginary part of the eigen- 
frequency implies an exponential growth in time of the 
corresponding excitation mode, leading to decay of the 
multiply-quantized vortex state. This spectral instability 
is a general parametric phenomenon occurring when sev- 
eral modes compete during coherent evolution and has 
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FIG. 3: Density dependence of the decay process. The time 
scale for the decay process of doubly quantized vortex states 
was measured by observing the vortex cores and classifying 
them as one vortex (open circles) or two vortices (solid cir- 
cles). Data were collected with three axial trap frequencies 
f z = 2.7, 3.7, 12.1 Hz and the interaction strength an z was 
controlled by changing the atom number by rf induced evap- 
oration before imprinting. Typical absorption images for (a) 
fast decay at low density (an z — 1.5) and (b) slow decay at 
high density (an z = 10.1). The field of view in the absorption 
images is 300 /im x 300 /im. (c) The separation of two visible 
cores vs. the hold time for 2 < an z < 3 (solid triangles) and 
6 < an z < 8 (open triangles). The solid and dashed lines in- 
dicate the diameter of one vortex core and of the condensate, 
respectively. 



been studied in many other nonlinear physical systems 
(see, e.g., Ref. 0,H2 and references therein). 

For a doubly-quantized vortex state in a cylindrically 
symmetric condensate, it was theoretically found that 
there are two excitation modes with a complex eigenfre- 
quency One of them is confined inside the doubly- 

quantized vortex core; the growth of this so-called "core" 
mode induces splitting of the original doubly-quantized 
vortex core into two separate singly-quantized vortex 
cores. The other mode, having the conjugate eigenfre- 
quency, grows with the core mode in order to conserve 
energy. In the low density limit, this mode corresponds 
to the co-rotating quadrupole mode, leading to oscilla- 
tions in the surface shape of condensates. We always 
observed that the surface of condensates changed into a 
quadrupole shape as the two cores appeared, as shown in 
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FIG. 4: Examples for the dynamic evolution after imprinting 
a doubly quantized vortex: (a) Surface Exicitation. Regular 
density modulation of the surface was observed after 51 ms 
hold time for an z = 1.8 (b) same as (a) with a contour line. 

(c) Crossing of vortex lines. 55 ms hold time and an z = 8.4. 

(d) same as (c) with guide lines for vortex lines. The field of 
view is 270 fim x 270 /xm. 



Fig.[2{a), and the ellipticity was larger at lower density. 

The dynamical instability of the doubly-quantized vor- 
tex state is related to the magnitude of the imaginary 
part of the complex eigenfrequency, and, according to 
the numeric calculation in Ref. nonvanishine; imag- 
inary part of the eigenfrequency appears at an z < 3 and 
an z ~ 12, showing a quasi-periodic behavior as a func- 
tion of the interaction strength, an z . The experiment 
showed a monotonic increase of the lifetime with no hint 
of periodic behavior. However, the calculated instability 
is not directly comparable to the observed lifetime. The 
imaginary part represents only the initial instability. Our 
criterion for decay was the observation of two separated 
vortex cores. It is possible that the dynamical instability 
changes after the doubly-quantized vortex state is signif- 
icantly perturbed 0, El- It would be helpful to have 
more inclusive calculations leading to a lifetime directly 
comparable with the experiments. 

What is the further evolution of the two cores? Some of 
the images at low density showed a regular surface mod- 
ulation, as in Fig. Ufa), which was not seen in clouds 
with a single core. This indicates that higher-order 
surface modes are excited during the coherent evolu- 
tion . However, their reproducibility was insufficient 
for a systematic study. Several images, especially those 
labelled as "undetermined", suggest that vortex lines 
crossed mm, as in Fig. Hfc). In °ur system, it was 
difficult to trace the positions of the two cores beyond 
80 ms hold time. 

In conclusion, we observed how a doubly-quantized 
vortex splits into a pair of singly-quantized vortices, and 
found higher stability at higher atom density. The topo- 
logical phase imprinting technique is unique in generating 
doubly- or quadruply-quantized vortex states jlfj 135^ ; a 
key feature is the rapid preparation of well-determined 
vortex states which gives access to their dynamical in- 
stabilities and coherent evolution. 
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